Azimuthal anisotropy in a jet absorption model with fluctuating initial geometry in 

heavy ion collisions 
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The azimuthal anisotropy due to path-length dependent jet energy loss is studied in a simple jet 
absorption model that include event by event fluctuating Glauber geometry. Significant anisotropy 
coefficients v n are observed for n — 1,2 and 3, but they are very small for n > 3. These coefficients 
are expected to result in a "ridge" for correlations between two independently produced jets. The 
correlations between the orientation of the n th -order anisotropy induced by jet absorption (<!>„) and 
the n th -order participant plane ("3?^) responsible for harmonic flow are studied. Tight correlations 
are observed for n = 2 in mid-central collisions, but they weaken significantly for n 7^ 2. The 
correlations are positive for n < 3, but become negative in central collisions for n > 3. The 
dispersion between $„ and $^ p is expect to break the factorization of the Fourier coefficients 
from two-particle correlation v„ tn into the single particle v n , and has important implications for the 
high-px ridge phenomena. 



PACS numbers: 25.75.Dw 



I. INTRODUCTION 



13 Recently, a lot ol attentions are focused on the study 

14 of the azimuthal anisotropy of the particle production 

15 in heavy ion collisions at the Relativistic Heavy Ion Col- 
i6 lider (RHIC) and the Large Hadron Collider (LHC). This 
17 anisotropy is usually expanded into a Fourier series: 
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with v n and <&„ represent the magnitude and direction of 
n th -order anisotropy, respectively. At low p^ , Vn is driven 
by the anisotropic pressure gradient associated with the 
initial spatial asymmetries, with more particles emitted 
in the direction of largest gradients [l| . Asymmetries giv- 
ing rise to non-zero v n are associated with either average 
shape (for n = 2) or shapes arising from spatial fluc- 
tuations of the participating nucleons 0t3- They can 
be characterized by a set of multi-pole components (also 
known as "eccentricities") at different angular scale, cal- 
culated from the participating nucleons at (r, (f>) [3j : 
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29 The orientations of the minor axis for each moment n, 

30 also known as the participant plane (PP) are given by 

P p atan2((r 2 sinn^>), (r 2 cosncj)}) it 
n n n 

31 When fluctuations are small and linearized hydrodynam- 

32 ics is applicable, each moment of the flow v n is expected 

33 to be independently driven by e„ along <f> pp = $„ 3]. 

34 This may not be true when the fluctuations are large, as 

35 the non-linear effects may lead to significant mixing be- 

36 tween harmonic flow of different order [a. In this paper, 



they are assumed to be the same to facilitate the study of 
the correlations between $„ of different physics origins. 

At high pt (pt > 10) GeV, the v n is understood to 
be driven by the path-length dependent energy loss of 
jets traversing the medium, with more particles emitted 
along the direction of shortest path-length (or direction of 
smallest jet attenuation) [3, Q . In contrast to flow which 
is sensitive to the global geometry manifested through the 
global evolution of the created matter, this anisotropy is 
more sensitive to the local density experienced by the jets 
as they traverse the matter. Nevertheless, since both flow 
and jet quenching are influenced by the same geometry, 
the directions of largest pressure gradient for flow and 
direction of smallest jet attenuation are strongly corre- 
lated. In fact, they are often implicitly assumed to be the 
same in many theoretical calculations 043 • An explicit 
study of the correlation between the two directions can 
help clarifying this assumption. 

In this paper, we estimate the high-px anisotropy co- 
efficients v n and their associated directions 4 r ^ p (QP 
stands for "quenching plane") using a simple jet absorp- 
tion model with event by event fluctuating Glauber ge- 
ometry. We study the correlations between ^ pp and 
\l/^ p , and discuss implications of these correlations for 
the interpretation of the "ridge" phenomena in two- 
particle correlations (2PC). 



II. MODEL 

We use a simple jet absorption/ Glauber model of 
Ref. [l5[ to calculate v n and V E , ^ P . This model has been 
used previously to study the centrality and path-length 
dependence of single particle suppression Raa, dihadron 
suppression /aa and V2- Back-to-back jet pairs are gen- 
erated according to the binary collision density profile 
(p c ) in the transverse (xy) plane with random orienta- 
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FIG. 1: (Color online) The complete set of output obtained in the jet absorption model for one event in 0-5% centrality interval: 
(a) the participant density profile (p v ) ; (b) the collision density profile (p c ), (c) the probability distribution of the path-length 
integral 7i, (d) the probability distribution of jet surviving the exponential attenuation, (e) the distribution of survival rate 
as function of azimuth angle, (f) the initial spatial asymmetry of the participants calculated via Eq. [7] The original impact 
parameter of the event is aligned along the i-axis. 



tion. They are then propagated through the medium, 
whose density is given by the participant density profile 
(p p ). Both profiles are generated with a Monte Carlo 
Glauber model with event by event fluctuation of po- 
sitions of nucleons in Au ions [16]. The nucleons are 
sampled from a Woods-Saxon distribution with a radius 
of 6.38 fm and diffuseness of 0.535 fm, with a nucleon- 
nucleon cross-section of a nn = 42 mb. In order to have 
smooth distributions for p c and p p , the nucleons are as- 
sumed to have a Gaussian profile in transverse plane with 
a width of ro = 0.4 fm in x and y direction similar to 
Ref [1Q|. The value of r$ is varied from 0.2-0.4 fm, and 
the nucleon is also assumed to be a uniform disk with a 
radius of -\/<7 nll /7r/2 = 0.58 fm. However the final results 
are found to be insensitive to the details of the nucleon 
shape, except in peripheral collisions. 

The jet quenching is implemented via exponential at- 
tenuation / = e _K/ , where the matter integral I is cal- 
culated as 
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for jet generated at r = (x,y) and propagated along 



91 direction v. They corresponds to l m+1 dependence of 

92 absorption (ex l m dl) in a longitudinal expanding or 1+1 D 

93 medium (oc l/(lo + 1)) with a thermalization time of Iq = 

94 cTrj. The Iq is fixed to by default, but we have checked 

95 the v n do not change much for Iq < 0.3 fm [lfj. The 

96 two cases, m = 1 and m = 2, are motivated for the I 

97 dependence expected for radiative and AdS/CFT energy 

98 loss in 1+1D medium [IT], [l8|], respectively. 

99 The absorption coefficient k controls the jet quenching 
wo strength and is the only parameter in this calculation. 
ioi It is tuned to reproduce Raa = (e~ Klm ) ~ 0.19 for 0- 

102 5% 7T° data at RHIC after averaging over many Glauber 

103 events [lj|. This leads to a value of ft = 0.1473 fm^ 1 and 

104 0.0968 fm -2 for m = 1 and 2, respectively. 



III. RESULTS 

Figure [1] summarize the basic information obtained 
from this procedure for one typical Au-Au event in 0-5% 
centrality interval. Panels (a) and (b) shows the den- 
sity profile for p p and p c , respectively. Panel (c) shows 
the normalized probability distribution of I±: P(Ii(<f))), 
which is obtained by calculating I\ over all possible di- 
jet production point p c and jet propagation direction (f>. 
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This distribution exhibit characteristic high density and 
low density regions in (Ii,(f>) space, presumably reflect- 
ing spatial correlation between the p c and p p profiles. 
Panel (d) shows the normalized probability distribution 



of the attenuation 



~ — kI\ 



Panel (e) shows the 



,— k!\\ 



eraged along the y-axis in Panel (d) as a function </>, which 
is precisely the azimuthal angle dependent suppression 
Raa (</>)• A clear anti-correlation can be seen between 
the peak magnitude of the RaaW and breadth of the 
I\ distribution in Panel (c). This distribution can also 
be obtained by randomly generating many di-jet pairs 
according the p c and propagating them through p p via 
Eq. [6] We expand it into a Fourier series: 



RaaW = #aa(1 +2V^ QP c, „ 



'E' 



(0-$Q p )), (6) 
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126 where R\a represents the average suppression, v® p and 

127 <f>Q p represent the magnitude and direction of n th -order 

128 harmonic of emission probability distribution, respec- 

129 tively. Similar studies of -Raa(^) were pursued before 
i3o in Ref. 12] for a pQCD energy loss in a event by event 

131 hydrodynamic underlying event. However it focused pri- 

132 marily on the influence of fluctuations on the event- 

133 averaged Raa(4>) distribution relative to the 2 nd -order 

134 event plane (EP). 

135 Figure Q] (f ) shows a distribution calculated from e„ 

136 and $S P : 
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137 It visualizes the shape of the initial geometry that is 
us transformed into the final momentum anisotropy via ei- 
139 ther flow or jet quenching. A good alignment is seen 
wo between $ pp and $^ p for n < 3. It also shows that the 

141 large e n for n > 3 are strongly damped after jet absorp- 

142 tion, leading to very small values of v® p for n > 3. 

143 The study shown in Fig. [T] can be repeated for many 

144 events. We divide the simulation data into 5% centrality 

145 intervals, each containing about 2500 events. Figure [5] 
we shows the distribution of $ pp — $^ p f or two centrality 
147 intervals. Strong positive correlations are obtained for 
us n = 1, 2 and 3, while the correlations are rather weak or 

149 even become negative for n > 3 1 . 

150 In heavy ion collisions at RHIC and LHC, the v n is 

151 usually measured from particle distribution relative to 

152 ty n via Eq. Q] [20| . However it has also been derived from 

153 the Fourier coefficients of two-particle correlation in rel- 

154 ative azimuthal angle A<fi = </> a — <fi h [21j : 
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with 



Vn,n(PT,Pr) = V n (p^)v n {p\) 



(8) 



(9) 



1 < t , JJ p in Eq. is calculated with r 2 weighting. We have also re- 
peated the study using r n weighting for n > 1 and r 3 weighing 
for n = 1 [jj, but very little differences are seen. 
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FIG. 2: (Color online) The correlation between participant 
plane $ pp and quenching plane $„ p for n = 1 — 6 calculated 
for 7i path-length dependence and for (a) 0-5% and (b) 20- 
25% centrality interval. 



156 The fact that the quenching plane and participant plane 

157 do not align exactly with each other implies that the 

158 v n measured relative to ^ pp is not the same as those 

159 contributing to the 2PC in Eq. [51 In other words, it is 
wo possible that the v n obtained from single particle analysis 
i6i is only a fraction of the true anisotropy resulting from jet 
162 quenching: 



= v QP 



<cosn(c£. pp -<i>Q p )} 



(10) 



Since what is measured in experiment is the event plane 
not the PP, it is important to check whether the event 
plane align with QP or not, for example in a hydrody- 
namic model calculation. 

Figure [3] (a) and (c) summarize the centrality depen- 
dence of v% p for n = 1-6 and for I\ and I2, respec- 
tively. Significant v® p signals are observed for n < 3, 
while higher-order v^ p are usually smaller than 1%. The 
V2 and 7J4 -Vq all show strong centrality dependence, 
while the ly and vjg show little centrality dependence 
100. Interestingly, the value of the v± is 



173 for N, 



part 



> 



,Q p 



174 consistently larger than that for 7J3 , and it even ex- 

175 ceeds iy value in most central collisions. This behavior 

176 suggest that the path-length dependence of energy loss 

177 and initial dipole asymmetry from fluctuations corrobo- 

178 rate to produce a large v^ . This large v± is expect 

179 to contribute to the high-pT «i signal observed by the 

180 ATLAS Collaboration [21(. Figure [3] also shows that the 
i8i I2 type of path-length dependence induces significantly 

182 larger v^ p than that for I\ : the increase is almost a fac- 

183 tor of two for n = 1 and n = 3. This is also observed in 

184 other studies before [TTj, [l8| . 

185 Figure [3] (b) and (d) summarize the centrality depen- 

186 dence of (cosn($ pp — $^ p )) for n = 1 — 6 and for I\ 
is? and Ii , respectively. As indicated by Eq. [lOl this repre- 

188 sents the reduction factor of the v n when it is measured 

189 relative to the < I >PP . The reduction is small for n = 2, 

190 except in central collisions where it reaches 15% for I\ 
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FIG. 3: (Color online) The centrality dependence of 
anisotropy coefficients t>„ (left panels) and correlation be- 
tween the participant plane and quenching plane (cos n($^ p — 
<E>^ P )) (right panels) for 7i type of path- length dependence 
(top panels) and I2 types of path-length dependence (bottom 
panels). Note that the values of v„ are positive by construc- 
tion according to Eq. [5] 



191 and 30% for 12- However the reduction is significantly 

192 larger for n = 1 and 3, reaching about 50% for n = 1 in 

193 mid-central collisions. The (cosn($^ p — $^ p )) value be- 

194 comes negative for n > 3 in central collisions, reflecting 

195 an anti-correlation between <I> PP and $^ p (already shown 

196 in Fig. [2]). Interestingly, (cosn($ pp — $^ p )) values for 

197 n = 1 are always smaller than that for n = 3 (more 

198 misalignment), while v^ is always larger than v$ . 

199 The dispersion between the $Q p and $ pp has impor- 

200 tant implications on the factorization relation Eq. [5] The 

201 factorization of v n>n into v n is obviously valid for cor- 

202 relations between two low pt particles (soft-soft corre- 

203 lation) as both are modulated around $ pp . The fac- 

204 torization should also be valid for correlation between 

205 a low-px particle and a high-px particle (soft-hard cor- 

206 relation) since it involves the projection of the v n onto 

207 $ pp , i.e. iv™(Pt>Pt) = v n (p^)v% p {p\){cosn(<S> pp - 
20a $^ P )) = Vnip^}Vn{p\) ■ Experimental data indeed sup- 

209 port this |2ll . |23| . However the correlation between two 

210 high-px particles from two independent hard-scattering 

211 processes (hard-hard correlation) is expected to be larger 

212 than the product of the two single particle v n : 



.(p*,^) = ^ p ( Pt )^ P (Pt) 



(11) 
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FIG. 4: (Color online) The expected long-range structures for 
correlations between two high pt particles from independent 
hard-scattering processes. They are shown for two typical 
events in 0-5% centrality interval (left panels) and 20-25% 
centrality interval (right panels); each of them should be re- 
garded as the distributions obtained for many events with 
identical initial geometry. 



_ Vn(PT) V n(j>T) 

<cosn($ pp -$Jr)> 2 ' 

213 Therefore, the factorization can not work simultaneously 

214 for soft-soft, soft-hard and hard-hard correlations. 

215 The large anisotropy coefficients v^ p also has im- 

216 portant consequences for the "ridge" observed in two- 

217 particle correlations |2lTJ23l |. This "ridge" is understood 

218 to be the result of the constructive contribution of har- 

219 monies at A</> ~ 0. In the literature, it is referred to 

220 as either the "soft-ridge" [2J, [25| for soft-soft correlation 

221 or "hard-ridge" |22l . [26j | for soft-hard correlation, respec- 

222 tively. Here we show that the correlation between two 

223 independently produced high-px jets can also produce 

224 the "ridge" -like structure. This "hard-hard ridge" can be 

225 calculated on a probability basis event- by- event by sim- 

226 ply self-convoluting the Raa {4>) distribution like Fig. Q] 

227 (e). Examples of these structures are shown in Fig. 0] 

228 for two representative events in both 0-5% and 20-25% 

229 centrality intervals. The magnitude of the ridge, as well 

230 as the away-side shape changes dramatically from event 

231 to event. They also changes a lot between the I\ and I2 

232 types of path-length dependence jet absorption. 

233 Figure [S] show the long-range structures (solid lines) 

234 obtained from the jet absorption model, averaged over 

235 many events. The ridge magnitude increases with cen- 

236 trality to about 1.5% (4%) for I\ (I2) path-length de- 

237 pendence in mid-central collisions. This signal should be 

238 measurable with the large statistics dataset from LHC. 

239 The dashed lines in Fig. [S] show the 2PC predicted from 

240 the v n measured relative to $ pp . Clearly the misalign- 

241 ments between $^ p and $ pp reduces the ridge magni- 

242 tude. The reduction is almost 50% in most central colli- 

243 sions, but decrease to about 20% in mid-central collisions. 




246 large and measurable. 



IV. CONCLUSION 

The anisotropy of high-px particle is studied in a sim- 
ple jet absorption framework with event by event fluctu- 
ating geometry. The harmonic coefficients v n are found 
to be significant for n = 1 — 3 (> 1%) but become very 
small for n > 3. The correlation between the quenching 
plane and participant plane are studied. A strong de- 
correlation is found for n = 1 in central collisions and 
for n — 1 and 3 over the full centrality range. The cor- 
relations become negative for n > 3 in central collisions. 
This de-correlation, if also confirmed between the event 
plane and the quenching plane (e.g via hydrodynamic 
model that has dijets embeded), is expected to break 
the global factorization of the two particle Fourier coef- 
ficient - 



FIG. 5: (Color online) The expected long-range structures for 
correlations between two high px particles from two indepen- 
dent hard-scattering processes (solid lines) and those calcu- 
lated from single particle i>„ relative to participant planes 
(dashed lines) for various centrality intervals. They are aver- 
age distribution over many events for a given centrality inter- 
vals. The thick (thin) lines denote the 7i (72) type of path- 
length dependence. 



2« This suggests the difference between the measured ridge 
245 and those predicted by the event plane method could be 
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into the v n for the two single particles. It 
would also imply that the high-px v n measured relative 
to the event plane could be significantly smaller than 
the true anisotropy from path-length dependent jet en- 
ergy loss. These jet quenching v n also give rise to long 
range "ridge" structure in two-particle correlations. The 
predicted ridge amplitude is on the order of 0.5-4% de- 
pending on the centrality and functional form of the I 
dependence of the energy loss, and should be measurable 
at the LHC using the correlations between two high-px 
particles with a large rapidity separation. Our study bear 
some similarities to Ref. [14| . However, Ref. [14j uses a 
cumulant expansion framework instead of Monte Carlo 
Glauber model for initial geometry, and that it focus on 
the soft-hard ridge instead of the hard-hard ridge in our 
case. 
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